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Modeling  and  damage  detection  for  complex  structures 


Challenges: 

•  Component-level  damage  affects  system-level  dynamics 

•  Fast  re-analysis  is  needed  to  reduce  computational  cost  of 
large-scale  finite  element  models 

•  Cracks  create  nonlinear  dynamics  (much  harder  to  tackle) 

•  Structural  health  monitoring 
information:  sensors 


_  ,  Vehicle  frame  model  (developed  by  Prof.  Hulbert,  Dr.  Ma,  Dr.  Hahn  of  the  Univ.  of  Michigan) 

Approach: 

•  Apply  component-based  methods  to  assemble  system-level 
reduced-order  models  (ROMs)  of  damaged  structures 

•  Employ  linear  approximations  of  nonlinear  (cracked)  structural 
dynamics 

•  Combine  above  into  sensor  placement  /  measurement  point 
selection  algorithm 
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Reduced  Order  Models:  Overview 


■  Dynamic  analysis  of  invariant  complex  structures 

>  Projection  by  lower  modes  of  the  large-scale  eigenvalue  problem 


■  Dynamic  analysis  of  damaged  complex  structures 

>  Projection  by  proper  basis  of  the  large-scale  eigenvalue  problem 

>  Proper  basis  can  be  defined  for  each  damage  type:  cracks,  dents  and 
other  structural  variations  of  complex  structures 
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Reduced  Order  Models:  Substructuring 


■  Assemble  ROMs  of  system  (e.g.,  frame)  from  finite  element 

analyses  of  components  and  subcomponents 


■  Efficiently  predict  vibration,  loading,  stress  in  critical  regions 


System  Level: 
Vehicle  Frame 


Finite  Element  Model 
of  Frame 


Dynamic  stress  for 
component  mode 
(left  rail) 


Component  Level: 
Left  Rail 


Subcomponent  Level: 

Rail  Sections, 
Reinforcement  Plates 


0,  Reinforcement 
^  Plates 
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Reduced  Order  Models:  CB-CMS 
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-  Superscript  C  :  Constraint  part 

-  Superscript  N  :  Internal  part 
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Reduced  Order  Models:  Parametric  Models  (PROMs) 


Enable  fast  re-analysis 

Subcomponent  dynamics  evaluated 
at  sampled  parameter  values 

System-level  response  expressed 
as  function  of  parameter  changes 


Global  PROM  (Parametric  Reduced  Order  Models) 

■  Balmes:  Collected  eigenvectors  at  sampled  points  in  the 
parameter  space 

Problem:  Overhead  computational  cost  to  get  the  modal 
matrix  to  project  the  FE  model 

-  CMB-PROM  (Component  Mode  Basis  PROM) 

■  Zhang  (2005):  Collect  fixed  interface  normal  modes  and 
global  interface  mode  and  project  the  FE  model. 

Problem:  Global  analysis  not  substructural  analysis 

-  Component  PROM 

Park  (2008):  Developed  PROM  for  substructural  analysis 
Problem:  a  single  design  component  is  tackled 


Nominal 


Reinforcement  plate 
thickness  changes 


Upper 

Bound 


Multi-component 
PROM  (MC-PROM) 
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Reduced  Order  Models:  Static  Mode  Compensation 


Geometrical  variations  of  the  structure  (dents) 

■  Lim  (2004):  used  SMC  for  vibration  of  turbomachinery 
bladed  disks  for  geometrical  mistuning  using  SMC 


/th  healthy 
system  mode 
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Effect  of  damage 
by  assuming  external  force 
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Global  structure  analysis 
^  not  component-level  analysis 
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Basis  shape  Damaged  blade 

mode  (34,563  Hz) 

Component  Mode  Synthesis 

with  Static  Mode  Compensation 

(SMC-CMS) 


Reduced  Order  Models:  Nonlinear  Dynamics:  Cracks 


Cracks  in  the  structure 

■  Crack  surfaces  open  and  close  during 
vibration:  nonlinear  vibration 

■  Hybrid  Frequency /Time  Domain 
method  (Poudou  2003) 

■  Bilinear  Frequency  Approximation 
(Shaw  1983):  no  mode  information 


V 
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Reduced  Order  Models:  Bilinear  Frequency  Approximation 


Exact  for  nonlinear  vibration  frequency  of  a  piecewise  linear  oscillator 


■  Bilinear  frequency  approximation  (BFA)  for  multiple  DOF  (Chati  et  al.,  1997) 

■  BFA  using  general  3D  finite  element  model  (Saito  et  al.,  2009) 
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Reduced  Order  Models:  Bilinear  Mode  Approximation 


Opening  mode  (I)  Sliding  mode  (ll)Tearing  mode  (III) 

Crack  modes 


■  Manage  boundary  conditions  on  the  crack:  open  and  closed  cases 

■  Crack  open:  open  boundary  condition:  DOF  on  crack  surface  are  free 

■  Crack  closed:  sliding  boundary  condition:  free  sliding  inside  crack  surface 

■  Mode  approximation:  shape  of  vibration  is  a  linear  combination  of  mode 
shapes  for  open  and  closed  crack  cases  (dominant  coherent  structures) 


Bilinear  Mode  Approximation  (BMA) 
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Reduced  Order  Models:  Framework 


Analysis  Framework 

Divide  the  global  structure  into 
substructures  with  or  without  damage 

■  Apply  Craig-Bampton  CMS  (CB-CMS)  for 
substructures  which  do  not  have  any 
damage  or  variability 

o 

■  Apply  MC-PROM  for  the  substructure  with 
model  variations  (e.g.  uncertainties) 

■  Apply  BFA  for  cracked  structure  analysis 


o  O 


Assemble  substructures  for 
M&S  of  system-level  response 
under  various  damage 
locations  and  crack  lengths, 
uncertainties, 
design  changes 


Core  technologies 

■  CB-CMS 

g  Multi-Component  PROM 

■  SMC-CMS 

■  Bilinear  Frequency  and  Mode  Approximation 


Efficient  framework  for 
damage  detection  and  for 
structural  predictions 
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Example:  L-Shape  Plate  :  Dents  and  Thickness  Variations 
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ROM 


Thickness  variations 


y  =  0.03  (structural  damping) 

u  :  physical  coordinates 
q  :  modal  coordinates 
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Substructure 

Thickness,  Case  1 

Thickness,  Case  2 

1 

0.4  mm  - 

->  0.473  mm 

0.4  mm  - 

->  0.435  mm 

6 

0.4  mm  - 

->  0.422  mm 

0.4  mm  - 

->  0.491  mm 

7 

0.4  mm  - 

0.493  mm 

0.4  mm  - 

->  0.481  mm 

i  cicaoc 


Velocity  (mm/sec)  Velocity  (mm/sec) 


Results:  L-Shape  Plate:  Forced  Response 


UNO 


Frequency  (Hz) 


Frequency  (Hz) 


Results:  L-Shape  Plate:  Dents,  Thickness  Variations  and  Crack 


Free  Response 


Ur 
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Results:  Vehicle  Frame  :  Dents  and  Thickness  Variations 


Each  reinforcement  frame  has  thickness  variation 
Engine  cradle  has  a  dent 


Substructure 

Thickness,  easel 

Thickness,  case2 

L_rear 

3.0378  mm  - 

4.6268  mm 

3.0378  mm  — > 

5.5788  mm 

L_front 
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Results:  Vehicle  Frame:  Forced  Response 
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Results:  Vehicle  Frame:  Dents,  Cracks  and  Thickness  Variations 


UNCLASSIFIED:  Distribution  A.  Approved  for  public  retease 


L  front 


Engine  Cradle 


L_rear 


Crack  length  varying 
along  the  surface 

88.89% 


11.11% 


Bilinear  frequency  (Hz)  Bilinear  frequency  (Hz) 


Results:  Vehicle  Frame:  Free  Response 
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Sensor  Placement  Algorithm  for  Cracked  Structure 


General  sensor  placement  algorithm:  EIDV 

-  Effective  independence  distribution  vector  (EIDV)  [Kammer,  1991 ;  Penny  et  al.,  1994] 

-  From  the  real  modal  matrix,  the  EIDV  algorithm  is  executed. 

Problem 

-  The  augmented  BL  modal  matrix  ^BLcan  be  linearly  dependent 

Solution 

-  Use  left  singular  vector  U  of  within  the  criteria  to  EIDV 

®Si:(MxN)  M  :  Number  of  candidate  measurement  DOF,  N  :  Number  of  mode 
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Algorithm  for  modified  EIDV  with  Left  Singular  Vector 


-  Calculate  the  mode  shape  for  the  healthy  and  damaged  structures 
for  open  and  closed  cases  in  reduced  order  domain. 

-  Construct  the  BL  modal  matrix  for  the  healthy  and  damaged 
structures. 

-  Find  the  sensitive  mode  shapes  (and  their  frequencies)  by  using  the 
generalized  MAC  matrix. 

-  Make  bilinear  augmented  modal  matrix  by  the  sensitive  mode 
from  the  modified  MAC  matrix. 

-  Obtain  the  left  singular  vector  U  of  <&BL  and  make  <bSVD  which  is 
consist  of  left  singular  from  ^  to  Up  based  on  the  criteria 
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^  ^  SVD  ^SVD' 


-  Calculate  Fisher  information  matrix  given  by 

-  Calculate  effective  independence  distribution  vector  (EIDV),  the 
diagonal  of 


E  =  <*>sro  A 


<I> 


SVD  ‘ 
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Example:  Cracked  plate 
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Mode  for  healthy  structure 
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1.  Calculate  the  mode  shapes  for  open 
and  closed  states  at  each  crack  length 

2.  Construct  the  BL  modal  matrix  for  the 
healthy  and  damaged  structures 

3.  Find  the  sensitive  mode  shapes  by 
using  the  generalized  MAC  matrix 


40%  cracked  structure 
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SVD  ratio 


SVD  ratio  versus  number  of  measurement  locations 


Assessing  EIDV  using  SVD  ratio 
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Results:  Measurement  point  selection 


Find  the  sensitive  mode  shapes  by  using 
the  generalized  MAC  matrix 


Apply  modified  EIDV  using  bilinear  modal 
matrix  assembled  with  sensitive  modes 
to  select  measurement  points 


(15  points  in  this  example) 


4th  to  9th  mode  shapes  are  sensitive  for  healthy  and  40%  cracked  structure 
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Summary 


Modeling  and  simulation  of  damaged  structures 

•  Reduced-order  models  for  dents,  thickness  changes,  etc. 

•  Fast  reanalysis  methods 

•  Bilinear  approximations  for  predicting  nonlinear  effects  of  cracks 

Sensor  placement  (measurement  point  selection)  method 

•  Bilinear  mode  approximation  (BMA) 

•  EIDV-based  algorithm  for  point  selection 

Future  work 

•  Applications  to  SHM  of  complex  structures,  joining/fastening 

•  Applications  to  design  for  reliability,  observability 
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